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Resumen 


LÓPEZ-CANTENS, J.G., PRADO-HERNÁNDEZ, V., DE 
LEÓN-MOJARRO, B. e RUIZ- CARMONA. V.M. Verificación 
de procedimientos numéricos de simulación de maniobras 


en compuertas. Tecnología y Ciencias del Agua. Vol. V, núm. 1, 
enero-febrero, 2014, pp. 43-56. 


La aplicación de las ecuaciones de Saint Venant en los 
tramos de canales de riego requiere definir condiciones de 
frontera interna en los sitios donde se ubican estructuras 
de control como compuertas, considerando su ecuación 
dinámica para transmitir información hidráulica hacia 
ambos lados de la estructura, por lo que en el trabajo se 
verifican tres procedimientos numéricos de ejecución de 
maniobra de compuertas, con la finalidad de representar el 
movimiento de las estructuras de control, minimizándose el 
desfasamiento de caudales posteriores a una ejecución. Para 
ello se utilizó el canal Rehbock de la Universidad Autónoma 
Chapingo, México, realizándose de forma automatizada la 
medición de los niveles de agua y gasto, obteniéndose que 
en la maniobra progresiva con reinicialización de tirante, 
la tendencia del gasto en la vecindad y la compuerta es 
similar, atenuándose las oscilaciones en el tiempo, el cual 
se corresponde con la dinámica del agua en los canales; 
además, cuando la variable a regular es el tirante, se 
recomienda usar la ejecución de maniobra instantánea 
o la progresiva con reinicialización de gasto; pero si la 
variable es el gasto, es preferible la maniobra progresiva 
con reinicialización de tirante; además, el procedimiento de 
simulación de la maniobra progresiva con reinicialización 
de tirante es el que mejor representa la dinámica del agua 
en la estructura. 


Palabras clave: movimiento de compuertas, modelo de 
simulación, régimen no permanente, estructuras de control. 


Abstract 


LÓPEZ-CANTEÑNS, J.G., PRADO-HERNÁNDEZ, V.,, DE 
LEÓN-MOJARRO, B. €; RUIZ- CARMONA. V.M. Verification 
of numeral procedures for simulation of gate maneuvers. Water 
Technology and Sciences (in Spanish). Vol. V, No. 1, January- 
February, 2014, pp. 43-56. 


Applying Saint Venant equations in reaches of irrigation canals 
requires defining internal boundary conditions at the sites where 
they are located control structures like gates, considering the 
dynamic equation to transmit hydraulic information from one 
side of the structure to continue on the other. In this study three 
numerical procedures for maneuvering gates are tested to verify 
whether they represent the movement of these control structures, 
minimizing flow phase shifts after a maneuver. To this end, the 
Rehbock canal of the Autonomous University of Chapingo was 
used. Automated measurements of water depth and flow were 
taken. Results show that in progressive maneuvering with 
restart of water depth, flow through the gate and near it tends to 
be similar, while oscillations diminish over time, a phenomenon 
which is characteristic of water dynamics in canals. Moreover, 
when depth is the variable to be regulated it is advisable to 
use instantaneous or progressive maneuvering with restart 
of flow. On the other hand, if the variable is flow, progressive 
maneuvering with restart of depth is preferred. Simulation of 
progressive maneuvering with restart of depth is the procedure 
that best represents water dynamics in the structure. 


Keywords: gates movement, simulation model, unsteady flow, 
control structures. 
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Introducción 


Los sistemas de conducción y distribución se 
integran con un conjunto de canales sobre los 
cuales se instalan las estructuras necesarias 
de regulación. El funcionamiento de estos 
sistemas se evalúa como eficiente si se lleva 
a cabo sin exceso y sin falta de agua, y sin 
alterar demasiado el caudal de las extracciones 
laterales. Para mantener los tirantes y 
caudales constantes dentro de un cierto rango 
es necesario ma-niobrar las estructuras de 
regulación. 

La configuración de los sistemas de canales 
es compleja, ya que en ellos se alojan estructuras 
de control estáticas (vertedores) y dinámicas 
(compuertas) de diferentes tipos; importantes 
estructuras de extracción, como las de los 
canales laterales, y una gran cantidad de tomas 
directas o tomas granja con pequeños caudales, 
así como estructuras de aportes laterales. La 
representación matemática de este conjunto 
de elementos se puede realizar mediante un 
modelo de flujo no permanente del tipo Saint 
Venant que no tiene solución analítica (De León 
et al., 2002). 

Las ecuaciones de flujo no permanente son 
válidas para estudiar el movimiento del agua 
en los canales de riego, pero es necesario definir 
condiciones de frontera, internas o externas, 
para representar las estructuras de control 
y medición, siendo necesario incorporar en 
el proceso de modelación las ecuaciones de 
éstas, teniendo en cuenta su funcionamiento 
(Malaterre y Baume, 1998). 

Manz (1998) menciona que con los modelos 
de simulación dinámica se pueden predecir las 
condiciones hidráulicas de flujo en cualquier 
punto dentro de los sistemas de canales 
abiertos, como una función de las características 
físicas y de operación, así como predecir y 
evaluar la secuencia y el tiempo de maniobra 
que deben realizarse al lograr las asignaciones. 
Asimismo, señala que estos modelos simulan 
dos tipos básicos de características hidráulicas: 
los componentes de los canales abiertos y las 
estructuras de control hidráulico. 


La mayoría de los modelos matemáticos 
simulan de forma adecuada el movimiento 
del agua en los canales debido a que satisfacen 
los criterios técnicos de conservación de la 
masa, estabilidad, precisión, convergencia 
1991; 


Sin 


y consistencia (Clemmens et al. 
Goussard, 2000; CEMAGREF 2000). 
embargo, presentan problemas para simular 
la propagación de gasto producida por la 
maniobra de estructuras de control, lo cual 
puede producir oscilaciones de tipo numérico 
que prolongan el tiempo de transición de 
un régimen permanente a otro, y además se 
refleja en la propagación desfasada de los 
gastos debido a la no correspondencia entre el 
comportamiento físico y la solución numérica 
del fenómeno. 

Las deficiencias anteriores hacen que sea 
necesario profundizar en el estudio de las 
variaciones temporales de las características 
hidráulicas inducidas por el movimiento de 
las estructuras de control. De León et al. (2007) 
presentaron tres procedimientos matemáticos 
que combinan las ecuaciones de Saint Venant 
en la vecindad de las compuertas con las 
ecuaciones de la estructura para transmitir 
información hidráulica a ambos lados de la 
misma, por lo que el trabajo tiene como objetivo 
verificar los procedimientos numéricos de 
simulación anteriores con la finalidad de 
minimizar el desfasamiento de caudales 
posteriores a una ejecución de maniobras en 
una compuerta radial. 


Representación matemática de la 
operación de canales de riego 


Cuando se varían las aberturas de las com- 
puertas o los caudales en las tomas laterales, 
la dinámica del agua en un tramo de canal de 
riego a superficie libre limitado por estructuras 
de control se puede representar mediante las 
ecuaciones de Saint Venant, que de acuerdo con 
Chow (1994) están integradas por la ecuación 
de continuidad: 
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RL Y (1) 
0X ot 


donde T es el ancho de la superficie libre del 
agua (m); y, el tirante del agua (m); Q, el gasto 
(més); q, el gasto lateral de aportación (+) o 
extracción (-), por unidad de longitud del cauce 
(més"m*) x, la coordenada en la dirección 
longitudinal del cauce (m); t, el tiempo (s). 

Y la ecuación dinámica, denominada tam- 
bién como ecuación de momento: 


9Q AQ), Y 
ES +34 =8A(8,-S;)+kqV (2) 


Donde S ¡es la pendiente de fricción (adim.); 
S,, la pendiente del fondo del canal (adim.); g, la 
aceleración gravitacional (ms?); V, la velocidad 
media del flujo (ms”); A, el área hidráulica (m2); 
k, un coeficiente que toma el valor de cero si q > 
0 y uno si q<0 (adim.) 

Las ecuaciones anteriores conforman un 
sistema hiperbólico de ecuaciones diferenciales 
parciales no lineales de primer orden en gasto 
y tirante. Es un sistema altamente no lineal 
que no tiene solución analítica y por lo tanto 
se requiere de procedimientos numéricos para 
su solución. 

Las ecuaciones (1) y (2) en los modelos 
de simulación por lo general se resuelven 
mediante un esquema implícito de diferencias 
finitas. En efecto, su solución discreta en cada 
sección del canal y a cada paso de tiempo es 
la ecuación recursiva lineal siguiente (Liggett y 
Cunge, 1975; CEMAGREE, 2000): 


R¡AQ, +S¡Ay; =1, (3) 

donde: 
AQ; = AAQ;¿1 +BAY ¡,1 +C (4) 
Ay; = DAQ;,1 +EAY ¡1 +F (5) 


R, S y T son coeficientes que varían en 
función de la condición límite; A, B, C, D, E 


y F, coeficientes que varían de un paso de 
tiempo a otro y en cada sección del dominio 
de integración; Ay y AQ, las incógnitas en 
las secciones ¡e i + 1; í, un índice que define 
la sección de cálculo dentro del espacio de 
discretización. 

Las ecuaciones (3), (4) y (5) relacionan las 
variables tirante y gasto en cada paso de tiempo 
para dos puntos consecutivos (i, ¡ + 1) definidos 
por la discretización del espacio. La ecuación 
recursiva (3) se puede resolver fácilmente 
mediante el método del doble barrido. El 
primer barrido se realiza de aguas arriba hacia 
aguas abajo para obtener los coeficientes R, S, 
y T, en el espacio de discretización mediante la 
ecuación tipo (3) de la forma siguiente: 


RiaAQua +S 1 AY jr = Tis (6) 
donde: 
Ri =R A, +S¡D; (7) 
Sia =R¡B,+S,E, (8) 
Ta "TRE -S, 09) 


Con el segundo barrido, que se realiza de 
aguas abajo hacia aguas arriba del tramo se 
calculan las variaciones de tirante y gasto en la 
sección ¡ mediante la solución de las ecuaciones 


(4) y (5). 
Condiciones de frontera 


Para aplicar cualquiera de los métodos de 
solución de las ecuaciones del régimen no 
permanente es necesario definir las condiciones 
iniciales y de frontera. La primera por lo general 
está representada por la curva de remanso en el 
instante inicial f = f, que se obtiene cuando el 
régimen es permanente gradualmente variado 
(Arteaga, 1997), mientras que la segunda se 
refiere a la definición del tirante, la velocidad o 
el gasto en los extremos superior e inferior del 
tramo, al considerar todos los tiempos desde el 
inicio del cálculo. 
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Liggett y Cunge (1975), y Husain et al. (1988) 
señalan que existen dos tipos de condiciones 
de frontera, las exteriores, ubicadas en los 
extremos de la longitud del canal considerado, 
y las interiores (CE), que se encuentran en las 
partes intermedias del mismo, ocasionadas 
por una discontinuidad geométrica o debido a 
la presencia de una estructura de control. De 
León et al. (2002) asevera que la formulación y 
el tratamiento de las condiciones de frontera 
) y al final (CF,,) del canal 
no presentan dificultades para su modelación. 
De León et al. (2007) y CEMAGREF (2000) 
señalan que la condición de frontera aguas 
arriba está representada por la demanda de 
riego, la cual es un hidrograma de entrada 


externa al inicio (CF 


arr: 


del tipo Q, = A(t), que se puede linealizar en la 
siguiente forma para solucionar el modelo de 
Saint Venant mediante la ecuación (3): 


Q,(t + At) - Q,(t) =f,(t + At) — f,(t) o también: 
AQ, =Aft) 


donde el subíndice O representa el inicio del 
canal. 

Con la condición de frontera al inicio del 
canal (aguas arriba) se inicia el primer barrido 
para calcular los coeficientes R, S, y T, ya que 
la ecuación (3) es recursiva. Así se tiene que los 
coeficientes R,, S, y T, al inicio del canal son: 


R,=1.0, 5, =0.0 y T,=Q(t + At)- Q(t) 


Con el primer barrido se llega al final del 
canal con un sistema de dos ecuaciones y dos 
incógnitas: 


ReAQ + Si AY y = Li (10) 


FAQs AY y) (11) 
donde m es el último tramo en que se discretizó 
el canal. 

La ecuación (11) representa la condición 
de frontera aguas abajo, que por lo general se 
representa mediante una relación conocida 


tirante-gasto, cuyos coeficientes R,, S, y T,, 
de la ecuación (3) adoptan la siguiente forma 


(CEMAGREEF, 2000; De León y Prado, 2006): 


m/ 


R,=10,5, =-¿ yT, -00 
A m 

Representación matemática de 

estructuras de control 


En los sitios donde se ubican las estructuras 
de control no son válidas las ecuaciones de 
Saint Venant, por lo que es necesario introducir 
dos ecuaciones para transmitir información 
hidráulica hacia aguas arriba (sección 2) y aguas 
abajo (sección 3); una asegura la continuidad 
del gasto y la otra describe la dinámica en la 
compuerta (figura 1). Por lo tanto, en estos 
puntos se debe resolver el siguiente sistema de 
ecuaciones (Cunge, 1988): 


R¿AQ, +5,Ay, =T, (12) 
AQ, = AQ,: ecuación de continuidad 
en la estructura (13) 
F(Q,,42,43,4)=0: ecuación dinámica 
en la estructura (14) 


donde a es la abertura de la compuerta 
(m); el subíndice 2 representa la sección de 
discretización aguas arriba de la compuerta y 
el 3 aguas abajo de la estructura. 

La ecuación (12) es el modelo de Saint 
Venant linealizado y representa la información 
hidráulica aguas arriba de la estructura. Del 
sistema de ecuaciones (12), (13) y (14) se obtiene 
la ecuación que representa la información aguas 
abajo de la estructura (De León et al., 2007): 

R¿AQ) +S3Ayz =T, (15) 

Los coeficientes R,, S, y T, se obtienen de la 

linealización de la ecuación (14): 


Ll0,v13,0)=0 (16) 
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Figura 1. Discretización de un tramo de canal de riego. 


Mediante un análisis en el lugar de las 
estructuras del modelo Saint Venant, De 
León (1990) obtuvo una ecuación diferencial 
ordinaria de segundo grado que permite 
estimar el comportamiento hidráulico en la 
vecindad de la estructura. Dicha ecuación se 
expresa de la siguiente forma: 


1(4QY 20/40 ¿LT 
Tldy) Aldy) 2? 


La solución principal de la ecuación (17) 


gA=0 (17) 


en dQ /dy se obtiene cuando hay movimiento 
de estructura y resulta igual a la ecuación 
siguiente: 


(18) 


donde c = (g4/T)'? es la celeridad de las ondas 
producidas por la maniobra de la estructura 
(ms”). 

En la ecuación (18), el signo positivo indica 
la relación gasto-tirante inmediatamente aguas 
abajo, mientras que el signo negativo se refiere 
a la relación aguas arriba. 

La ecuación (17) sirve como base para 
establecer de 


procedimientos simulación 


de ejecución de maniobras que permitan 
minimizar el error de propagación de gastos. 
De esta manera, De León et al. (2007) plantea 
que en los puntos donde se localiza una 
singularidad se puede aplicar la ecuación (18) 
que identifica a la ecuación recursiva de tipo (3) 
si R,=1.0, S,=-(v + c)T y T,=0.0. Esta ecuación 
representa las líneas características que parten 
con información hidráulica hacia aguas arriba 
y aguas abajo de la estructura, e interceptan la 
malla de discretización de diferencias finitas. 


Ejecución de la maniobra de compuertas 
de una red de canales de riego 


La maniobra de compuertas de regulación 
en una red de canales de riego es una acción 
necesaria para mantener el gasto y el tirante 
cercano a valores de referencia previamente 
definidos (Litrico y Fromion, 2006). 

Si la estructura es estática, una aproxi- 
mación tangencial de la ecuación dinámica de 
la compuerta puede utilizarse para continuar 
la solución aguas abajo de la misma. Sin 
embargo, De León y Exebio (1990) señalan 
que cuando se introducen variaciones tem- 
porales a la abertura de las compuertas, el 
procedimiento de aproximación tangencial 
conduce a errores, los cuales son acumulativos, 
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puesto que en cada paso de tiempo, o a cada 
instante de variación, el error de ajuste entre la 
aproximación tangencial y la función exacta de 
la singularidad se agrega al nuevo valor de las 
variables de cálculo (gasto y tirante). 

De León (1990) señala que el problema de la 
aproximación tangencial requiere una forma de 
previsión de la variación de los niveles aguas 
arriba (dy) y aguas abajo (dy, ,) ante un cambio 
de la demanda de riego en la estructura (40). 
Lo anterior se logra mediante la ecuación (18), 
aplicándose el procedimiento siguiente: 


_ AQ 
dy; 3 (o, -c,)T, (19) 
ja HB 20 
Via (Dj1 +0 54) La | 


y así se obtienen los nuevos valores de los 
tirantes: ye =y; +dy; y ye = Ya + AY y donde 
n representa la etapa de tiempo previa a la 
variación del gasto. 

La ecuación recursiva tipo (3) representa un 
conjunto de rectas que describen la variación 
temporal del tirante (dy) y el gasto (4Q) en 
cada malla de la discretización del espacio. 


Al superponer en un plano (dy, dQ) las rectas 


anteriores y las curvas de abertura de compuerta 
como se muestra en la figura 2, permitieron a 
De León et al. (2007) analizar las trayectorias 
descritas por distintos procedimientos de 
ejecución de maniobras. 

Como en la ecuación del tipo (3), la variación 
temporal del gasto representa el aporte y la 
variación del tirante el almacenamiento, se 
debe cumplir para cada paso de tiempo la 
siguiente relación: aporte + almacenamiento = 
constante. 

Según la relación anterior al ejecutar 
una maniobra existen tres posibilidades: a) 
modificar el almacenamiento y el aporte; b) 
modificar el almacenamiento, y c) modificar el 
aporte. 

En la figura 2 se muestran las trayectorias 
descritas por los procedimientos de simulación 
de maniobras: el recorrido 13 corresponde a la 
posibilidad “a”; el recorrido 143 corresponde 
a la posibilidad “b”; y el recorrido 123 
corresponde a la posibilidad “c”. 

Según De León y Prado (2006), el tratamiento 
numérico de estas tres posibilidades genera 
dos procedimientos generales de simulación: 
maniobra instantánea obtenida con la posibi- 
lidad “a” y maniobra progresiva obtenida con 


mM 2” 


las posibilidades “b” y *c”. 


n+1 n+1 pa 2 
Ay y “Y 


R 


dle ro tds + +1, 1+1_ pan 


S 
(Después de la maniobra) 


NA n MA n n 
RQ +S%y=T 
(Antes de la maniobra) 


Figura 2. Representación de la maniobra de una estructura hidráulica. 
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La maniobra progresiva se denomina así 
porque al suceder un cambio en la demanda de 
riego se llega al punto de referencia (número 
3) durante un paso de tiempo posterior a 
la maniobra. Mientras que en la maniobra 
instantánea el desplazamiento al punto de 
referencia se realiza en el mismo instante que 
se ejecuta la maniobra. 


Materiales y métodos 
Caracterización de las instalaciones 


Para el desarrollo de la fase experimental se 
utilizó el canal Rehbock ubicado en el Labo- 
ratorio de Hidráulica del Departamento de 
Irrigación en la Universidad Autónoma Cha- 
pingo (UACh), en Chapingo, Estado de México. 

El canal Rehbock utilizado es de sección 
rectangular, con estructura metálica y paredes 
de cristal, tiene una longitud de 20 m, plantilla 
de 0.6 m, altura de 1 m y pendiente de 0.0005. Se 
alimenta a través de una válvula de 12” (0.305 
m), que entrega a un estanque amortiguador 
para evitar turbulencia a la entrada del canal, 
limitado al inicio por un vertedor de pared 
delgada sin contracciones laterales y con una 
altura suficiente para que no trabaje ahogado y 
una compuerta deslizante. 

Al final del canal se encuentra una com- 
puerta deslizante que trabaja como vertedor 
rectangular movible, y que permite controlar 
el nivel del agua en el canal. A 5.5 m del inicio 
se colocó y operó una compuerta radial de 
0.75 m de altura, 0.60 m de ancho y radio de 
compuerta de 0.9 m. 

La medición de los niveles de agua se realizó 
con un sistema automatizado, que consta de 
seis sensores continuos, cinco ubicados en toda 
la longitud del canal y uno aguas arriba del 
vertedor para determinar la carga y con ésta 
calcular el gasto de entrada. 


Procedimientos de simulación de maniobras 
de compuerta 


La maniobra de una compuerta en un canal de 
riego provoca el cambio de un gasto inicial Q, a 


un gasto final Q, . Esta situación se establece con 
una variación de la abertura de la compuerta 
de un valor inicial 4, a un valor final a, . 

Tomando en cuenta lo anterior, De León y 
Prado (2006) y De León et al. (2007) propusieron 
una metodología para simular maniobras de 
las estructuras de regulación que representen 
el fenómeno de propagación de gastos pro- 
vocados por una variación de la demanda. 
Estas metodologías consisten en reinicializar 
los valores del gasto, abertura y tirantes en la 
compuerta. La reinicialización de los valores 
de tirantes sólo se efectúa inmediatamente 
aguas arriba y aguas abajo de la compuerta; 
en el resto de las secciones, el modelo 
general propaga el efecto de la maniobra. La 
metodología propuesta por estos autores se 
detalla a continuación. 


Maniobra progresiva con reinicialización del gasto 


Este procedimiento consiste en reinicializar el 
gasto en la estructura, el cual corresponde a la 
variación de la demanda, así como al valor de la 
abertura, manteniéndose constantes los valores 
de los tirantes aguas arriba y aguas abajo. 

Se parte de que en un instante t = n se co- 
noce el estado hidráulico de la compuerta 
(Ay, AQ," a) y la ecuación de la recta tipo (3) 
que pasa por ese punto. Luego, mediante la 
maniobra, se puede instalar un gasto Q, = Q, 
+ AQ, para lo cual se requiere pasar del punto 
(Ago, AQ% 4) al puto (y, AQ.) 
(figura 2). 

Los autores señalan que si sólo se reinicializa 
el valor del gasto y los tirantes permanecen 
constantes, en realidad se pasa del punto (Ay, 
AQ”) al punto (Ay, AQ,'**), cuya ordenada es 
superior al valor deseado (AQ,”*). 

El punto (Ay, AQ,”**) corresponde a la 
intersección del sistema formado por las 
ecuaciones siguientes: 


Ay? = constante (21) 


AQ = f( Aya") (22) 
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Como la intersección de este sistema de 
ecuaciones no corresponde al punto de interés 
(Ay,"**, AQ”), los autores señalan que con el 
reinicio del gasto en un valor AQ,”** inferior 
al AQ,”*! se obliga al modelo a descender del 
punto (Ay, AQ,”*) hacia el punto de interés 
(Ay,"**, AQ,”*!) sobre la curva de la compuerta 
ara, 

Con este procedimiento se desplaza la 
recta del tipo (3) del punto (Ay, AQ”) al 
punto (Ay,”**, AQ.”**) en un paso de tiempo. La 
evaluación de los coeficientes R, S y T permite 
transmitir información hacia aguas abajo sobre 
la magnitud de la maniobra. 


Maniobra progresiva con reinicialización del tirante 


Este procedimiento consiste en reinicializar 
los valores de los tirantes aguas arriba y 
aguas abajo, y la abertura de la compuerta, 
manteniéndose constante el valor del gasto en 
la estructura. Las variaciones de los tirantes se 
estiman mediante las ecuaciones (19) y (20). 

Al igual que en el procedimiento anterior, 
se requiere pasar del punto (Ay”, AQ”, a” 
al punto (Ay,”**, AQ”, a,”*) (figura 2). Los 
autores señalan que si se reinicializan los 
valores de los tirantes aguas arriba y aguas 
abajo de la compuerta y el gasto permanece 
constante, se pasa del punto (AyAQ”) al 
punto (Ay,”*!, AQ”), cuya abscisa es inferior al 
valor reinicializado (Ay,”**). 

El punto (Ay,”*, AQ”) corresponde a la 
intersección del sistema formado por las 
ecuaciones: 


AQ” = constante (23) 


AQ; = f(ay2* aj") (24) 

La intersección de este sistema da como 
resultado un gasto inferior al demandado 
sobre la curva de la compuerta a,”*; luego, 
como el valor del tirante se reinicializa en un 
valor Ay,”**' superior al valor Ay,”**, se establece 
el sentido de desplazamiento de las variables 


del punto (Ay,”*, AQ”) hacia el punto de 
interés (Ay,”**, AQ.”**) sobre la misma curva de 
la compuerta a,”*!. 

Con este procedimiento, la recta del tipo 
(3) se desplaza de manera progresiva durante 
un paso de tiempo (Af) hacia el punto (Ay,”**, 
AQ,''"). De igual manera, el gasto pasa 
uniformemente del valor AQ” al valor AQ”** 
durante un paso de tiempo. 


Maniobra instantánea 


Este procedimiento consiste en reinicializar 
simultáneamente los valores del gasto y los 
tirantes en la estructura, obtenidos con las 
ecuaciones (19) y (20); así como la abertura 
de la compuerta. Los autores señalan que este 
procedimiento equivale a desplazar la recta 
del tipo (3) de manera instantánea del punto 
(Ay, AQ”) hacia el punto de interés (Ay,”", 
AQ,”**) (figura 2). Este punto corresponde a la 
intersección del sistema de ecuaciones: 


progre ¿Ae eE qu (25) 


aqu a Fay" ap) (26) 
donde los coeficientes de la ecuación (25) se 
determinan con la ecuación (18). 

En el proceso de simulación se consideró 
un intervalo de distancia Ax = 0.5 m y un 
intervalo de tiempo At = 10 s, que garantizan 
una representación adecuada de la dinámica 
del sistema. Además, teniendo en cuenta los 
criterios sobre el factor de ponderación del 
tiempo discutidos por Liggett y Cunge (1975), 
Saavedra y López (1996), y CEMAGREF (2000), 
se seleccionó un valor de 0.66 para garantizar 
la convergencia, precisión y estabilidad del 
modelo. 


Resultados y discusión 


Verificación de los procedimientos numéricos 
de simulación de maniobras en compuertas 


Para verificar los procedimientos numéricos de 
ejecución de maniobras en estructuras de control 
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propuestos por De León et al. (2007) se usaron 
datos experimentales del canal Rehbock y un 
modelo de simulación de flujo no permanente 
desarrollado por Prado (2006). Primero se 
calibró el modelo en régimen permanente para 
un gasto de operación de 0.095 més”, con el fin 
de obtener las condiciones iniciales. Durante el 
proceso se ajustaron los valores del coeficiente 
de rugosidad en el canal y el coeficiente de 
gasto en la estructura de control, siendo éstos 
de 0.009 y 0.6, respectivamente. 

Al analizar el perfil de la superficie libre 
del agua, en el cuadro 1 se observó que los 
residuales de los tirantes medidos y simulados 


a lo largo del canal son inferiores a + 0.004 m, 
mientras que el mayor error porcentual fue de 
0.67%, el cual se encuentra dentro de los rangos 
de error que consideraron como aceptables 
Arteaga (1997) y CEMAGREF (2000), que 
están alrededor del 1%. Por lo tanto, se puede 
considerar que el modelo es capaz de simular 
con suficiente precisión el comportamiento del 
agua en el canal para la condición analizada. 
Con la finalidad de observar el compor- 
tamiento de la ecuación tipo (3) se realizó 
una variación de gasto a la entrada del canal 
de 0.023 m's* sin efectuarse maniobra en la 
estructura y resultó la familia de rectas paralelas 


Cuadro 1. Calibración del modelo de simulación en el canal Rehbock. Régimen permanente. 


Distancia (m) n E 
Medido Simulado Residual 

0 0.009 0.689 

25 0.009 0.688 0.691 -0.003 
45 0.009 0.694 0.691 0.003 
5D 0.009 0.692 

9 0.009 0.595 0.595 0 
Sto 0.009 0.595 0.597 -0.002 
18 0.009 0.595 0.599 -0.004 
20 0.009 0.600 

0.025 


0.020 


0.015 


AQ (més?) 


0.010 


0.005 


0 0.02 


AQ5%, +0.3392Ay 2% =0.023 


AQ +0.3392 Ay ¿ ; =0.005 


0.04 0.06 0.08 


Ay (m) 


Figura 3. Familia de rectas aporte-almacenamiento en la estructura. 
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mostrada en la figura 3. Esta familia de rectas 
representa la variación del gasto y el tirante en 
la estructura hasta alcanzar un nuevo régimen 
permanente en 27 pasos de tiempo (4.5 min). 

Las rectas resultantes de la ecuación 
tipo (3) en la vecindad de la estructura, 
sobrepuestas a las curvas carga-gasto para 
diferentes aberturas, permitieron analizar los 
procedimientos de simulación de ejecución de 
maniobras propuestos por De León et al. (2007) 
para satisfacer una variación en la demanda de 
riego. 

Para satisfacer la variación de la demanda 
de riego se deben maniobrar las estructuras, lo 
cual trae consigo una variación de la geometría 
y de las características de flujo que dificultan 
el análisis de los términos que conforman las 
ecuaciones de Saint Venant, por lo que se realizó 
un análisis global del fenómeno posterior a una 
maniobra en una vecindad muy próxima a la 
compuerta. 

Para analizar el fenómeno, considerando 
descarga ahogada, las características geomé- 
tricas e hidráulicas en la estructura y su 
vecindad antes de la maniobra (t = 1) fueron: 


Abertura =0.19m  Tirante aguas arriba = 0.694 m 
Gasto = 0.095 m's*  Tirante aguas abajo = 0.595 m 
Para satisfacer una variación en la demanda 
de 0.023 m*s* un paso de tiempo después de la 
maniobra (t=1 +1) se obtuvieron las siguientes 
condiciones geométricas e hidráulicas: 


Abertura = 0.28 m 
Tirante aguas 
arriba = 0.667 m 


Tirante aguas abajo = 0.600 m 
Tipo de descarga = ahogada 


La abertura de la compuerta se obtuvo 
mediante su ecuación dinámica, mientras 
que los tirantes aguas arriba y aguas abajo se 
calcularon con las ecuaciones (19) y (20). 

Al superponer la familia de rectas de aporte 
—almacenamiento con las curvas de carga— 
gasto correspondiente a la estructura, se obtu- 
vo la figura 4, que representa las trayectorias 
que siguen los métodos de ejecución de 
maniobras de compuerta propuestos por De 
León et al. (2007). Del análisis de la figura 
anterior se observa que cuando se reinicializa 
el gasto y el tirante permanece constante (a), 


0.20 


0.16 


AQ (més) 


0.08 


0 0.05 


Después de la maniobra 
AQ 0.36574y "9! 0.1425 


Antes de la maniobra 
AQ” +0.3657Ay" =0.1307 i 


0.10 0.15 0.20 


Carga en la compuerta Ay (m) 


Figura 4. Recorrido del gasto y el tirante para la maniobra de la compuerta correspondiente a la variación de demanda. 
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la intersección de las ecuaciones (21) y (22) da 
como resultado un gasto de 0.141 m's* sobre 
la curva de la compuerta a” = 0.28 m, que 
es superior al de referencia (0.118 mf*s”). Sin 
embargo, con la reinicialización del gasto se 
impone el valor demandado, obligándose al 
modelo a descender hasta el punto deseado, 
cuyas coordenadas son Ay = 0.067 m y AQ 
= 0.118 m*s1, sobre la misma curva de la 
compuerta en un paso de tiempo. 

Para el caso en que se reinicialice el tirante 
y el gasto permanezca constante (b), se tiene 
un punto de intersección de las ecuaciones 
(23) y (24) que representa una carga en la 
compuerta de 0.044 m, la cual es inferior al 
valor de referencia (0.067 m) sobre la curva 
de la compuerta a;*' = 0.28 m. Esto ocasiona 
un gasto de 0.095 m's*, inferior al gasto 
demandado de 0.118 més”. La diferencia en el 
gasto se corrige con la reinicialización de los 
tirantes, imponiendo la carga necesaria para 
lograr el punto de referencia (0.067, 0.118) 
sobre la misma curva de la compuerta en un 
paso de tiempo. 

Cuando se realiza la maniobra instantánea 
de la compuerta ocurre una reinicialización del 
gasto y el tirante (c), por lo que el modelo se 
desplaza instantáneamente del punto inicial 
(0.097, 0.095) en la curva de la compuerta 


aí = 0.19 m hasta la ubicación del punto de 
referencia (0.067, 0.118) sobre la curva de la 
compuerta ay*' = 0.28 m. 


Evolución del caudal en maniobra progresiva 
Reinicialización del tirante 


Al imponer el tirante en el proceso de cálculo 
durante la simulación de la maniobra de la 
compuerta, se observó en la figura 5 que a 
partir de += 1 min, momento en que se efectúa 
la variación de la demanda el gasto aguas 
arriba, aguas abajo y en la compuerta tienen la 
misma tendencia, es decir, las oscilaciones se 
propagan al mismo tiempo y se amortiguan de 
la misma forma. Sin embargo, el modelo estima 
un paso de tiempo después de la variación de 
gasto, siendo ésta de 0.011 més”, valor que 
es inferior al 0.023 més* correspondiente al 
cambio de la demanda, por lo que se produce 
una subestimación del aporte ocasionado 
por la variación de almacenamiento; com- 
portamientos similares fueron encontrados 
por De León y Prado (2006), y De León et al. 
(2007). 

Las del gasto 
provocadas por la maniobra de la compuerta 


oscilaciones numéricas 


0.110 


0.105 


0.100 


Q (més) 


0.095 


0.090 


0.085 
0.0 0.3 0.6 0.9 1.2 L5 
Tiempo (min) 


En la compuerta 


0.5 m aguas abajo de la 
compuerta 


0.5 m aguas arriba de la 
compuerta 


21 24 


Figura 5. Maniobra progresiva con reinicialización de tirante y gasto constante. 
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Figura 6. Maniobra progresiva con reinicialización de gasto y tirante constante. 


se atenúan de forma rápida, ya que el tiempo 
en alcanzar un nuevo régimen permanente es 
de 1.8 min. 


Reinicialización del gasto 


Cuando se reinicializa el gasto en el proceso de 
simulación, se aprecia en la figura 6 que a partir 
de t = 1 min, momento en que se efectúa el 
cambio en la demanda, hay un desfasamiento 
del gasto en la estructura con respecto a sus 
vecindades aguas arriba y aguas abajo. La 
maniobra repercute de manera instantánea en 
la variación del gasto en la compuerta, mientras 
que en su vecindad sucede el cambio un paso 
de tiempo después, propagándose en fase las 
oscilaciones en ambos lados de la estructura. 
Además existe una diferencia de gasto de 0.011 
m's* entre la compuerta y las secciones aguas 
arriba y aguas abajo, que se atenúa a los 1.8 
min. 

La variación del comportamiento del gasto 
en el tiempo para las secciones analizadas se 
debe a la incompatibilidad entre la escala 
numérica y la escala física del fenómeno de 
propagación, coincidiendo con De León et al. 


(2007), lo cual pudiera solucionarse reduciendo 
el paso de tiempo en el proceso de simulación, 
pero traería como desventaja el incremento del 
tiempo de procesamiento en la computadora. 


Evolución del caudal en maniobra 
instantánea 


Cuando se realiza la reinicialización del gasto y 
el tirante, se observa en la figura 7 que a partir 
de que comienza el régimen no permanente (t = 
1 min) existe un comportamiento similar de los 
gastos con respecto a la maniobra progresiva 
que reinicializa el gasto, pero la amplitud de las 
oscilaciones es mayor. Como puede apreciarse, 
existe variación del gasto en la estructura al 
momento de realizar la maniobra, mientras 
que aguas arriba y aguas abajo se observa un 
paso de tiempo después y se propaga en fase, 
estabilizándose las oscilaciones numéricas en 
1.8 min. 

Como se reinicializa el valor del gasto, 
es de esperarse que al momento de efectuar 
la maniobra, el gasto en la estructura sea 
igual al demandado, como se observa en la 
figura 7. También se aprecia que es el único 
procedimiento de ejecución de maniobras en el 
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Figura 7. Maniobra instantánea con reinicialización de gasto y tirante. 


que desaparecen las oscilaciones numéricas en 
la compuerta. Sin embargo, aunque el tiempo 
en que ocurre el fenómeno transitorio es el 
mismo en los tres procedimientos, en éste las 
oscilaciones numéricas son mayores. 

Se verificaron los tres procedimientos 
de simulación de maniobras de compuerta 
propuestos por De León et al. (2007), los cuales 
minimizan el problema del desfasamiento 
entre el gasto de las compuertas consecutivo 
a la operación y los gastos en la vecindad 
de ésta, además pueden integrarse estos 
procedimientos en la ejecución de maniobras 
de estructuras de control en modelos de flujo 
no permanente. Sin embargo, la selección del 
método de operación dependerá del problema 
a resolver; en caso de que la variable a regular 
sea el gasto, entonces se debe usar la maniobra 
progresiva con reinicialización de los tirantes, 
pero cuando la variable a regular es el tirante, 
se debe utilizar la maniobra progresiva con 
reinicialización de gasto. 

Al evaluar el comportamiento de los 
procedimientos de simulación de operación de 
compuerta (figuras 5, 6 y 7) se pudo determinar 
que la maniobra progresiva con reinicialización 
de tirante es la que más se corresponde con el 
funcionamiento real de una compuerta en un 


sistema de riego, ya que las variaciones del 
gasto ocurren al mismo tiempo en la estructura 
y en su vecindad al momento de efectuar una 
maniobra. 


Conclusiones 


Se cumplieron las hipótesis planteadas por De 
León et al. (2007) sobre la aplicación de tres 
procedimientos de simulación de maniobras 
de compuerta en canales de riego. 

Los procedimientos de simulación de 
maniobras evaluados minimizan el problema 
del desfasamiento de gastos en la vecindad de 
las estructuras provocados por una variación 
en la demanda. 

El procedimiento de simulación con reini- 
cialización de tirante es el que mejor representa 
la dinámica del flujo ante una ejecución de 
maniobra de la compuerta. 
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